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Geoscientific 
Model Development 


Abstract. The Palaeoclimate Modelling Intercomparison 
Project has expanded to include a model intercomparison for 
the mid-Pliocene warm period (3.29 to 2.97 million yr ago). 
This project is referred to as PlioMIP (the Pliocene Model In- 
tercomparison Project). Two experiments have been agreed 
upon and together compose the initial phase of PlioMIP. The 
first (Experiment 1) is being performed with atmosphere- 
only climate models. The second (Experiment 2) utilises 
fully coupled ocean-atmosphere climate models. Follow- 
ing on from the publication of the experimental design 
and boundary conditions for Experiment 1 in Geoscientific 
Model Development, this paper provides the necessary de- 
scription of differences and/or additions to the experimental 
design for Experiment 2. 


1 Introduction 

The mid-Pliocene warm period (mPWP) is defined by the 
United States Geological Survey’s PRISM Group (Pliocene 
Research Interpretation and Synoptic Mapping; http:// 
geology.er.usgs.gov/eespteam/prism/index.html) as the inter- 
val between 3.29 and 2.97 Ma (according to the geomagnetic 
polarity timescale of Berggren et al. 1995), lying between 
the transition of oxygen isotope stages M2/M1 and G19/G18 
(Shackleton et al., 1995), in the middle part of the Gauss Nor- 
mal Polarity Chron (Dowsett et al., 1999). The “Time Slab” 
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represents a climatically distinct period during the Pliocene 
when Earth’s climate was, on the whole, warmer than present 
(Dowsett et al., 1999; Dowsett, 2007). 

The mPWP has been the subject of intense study for the 
last two decades. There are many reasons for this, but an im- 
portant driver has been our desire to understand the dynamics 
of past warm climates as a potential guide to understanding 
climate change in the future (Haywood et al., 2009). The 
mPWP is well suited to this task. The climatic signal (change 
from modern) is sufficiently large, for many geographical re- 
gions, to be differentiated from the noise generated by the 
uncertainties and limitations inherent in the techniques used 
for palaeoclimatic/palaeoenvironmental reconstruction. The 
interval was the last time in Earth history when global tem- 
peratures were significantly warmer than modern, over a pe- 
riod longer than any Quaternary interglacial. It is unique 
in that continental configurations were relatively unchanged 
from today, and geological proxies are superior to those of 
preceding warm periods due to improved geographic cover- 
age, more reliable biota-environment correlations and higher 
resolution stratigraphy (Dowsett, 2007). 

For the initial phase of the Pliocene Model Intercompari- 
son Project (hereafter referred to as PlioMIP), which is a sub- 
component of the wider/larger Palaeoclimate Modelling In- 
tercomparison Project (PMIP; Braconnot et al., 2007a,b) two 
experiments were agreed upon. The first is an experiment us- 
ing atmosphere-only climate models (hereafter referred to as 
Experiment 1), the details of which can be found in Haywood 
et al. (2010; this journal and issue), whilst the second experi- 
ment (hereafter referred to as Experiment 2) will utilise cou- 
pled ocean-atmosphere climate models. Both experiments 
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use versions of the US Geological Survey’s PRISM Group 
boundary condition data sets (e.g. Dowsett, 2007). This Spe- 
cial Issue of Geoscientific Model Development represents the 
first set of co-ordinated publications from PlioMIP. It (a) de- 
scribes the chosen experimental design for Experiments 1 
and 2, (b) includes a detailed description of the boundary 
conditions used in both experiments, and (c) presents con- 
tributions from each participating model group, describing 
how the boundary conditions were implemented into the dif- 
ferent climate models and the basic results from the exper- 
iments themselves. This detailed record for the rationale 
and specifics of the experimental design, construction of the 
boundary conditions data sets, and critically, how these were 
implemented into each climate model, will provide an invalu- 
able reference when the intercomparison phase of PlioMIP is 
reached. This will help the PlioMIP/PMIP community to un- 
derstand more easily the differences which will inevitably be 
observed between mid-Pliocene warm period (mPWP) sim- 
ulations. The purpose of this paper is to briefly describe the 
experimental design and boundary conditions for PlioMIP 
Experiment 2, focussing only on the aspects that differ from 
Experiment 1 . Participating groups should refer to Haywood 
et al. (2010), for full details of Experiment 1 and for details 
relevant to Experiment 2 that are not repeated here. 

2 Experimental Design - Experiment 2 

2.1 Integration, atmospheric gases/aerosols, solar 
constant/orbital configuration 

The experimental design for Experiment 2 is summarised in 
Table 1 . The experiment integration length was set to 500 yr 
in accordance with CMIP5 guidelines (Coupled Model Inter- 
comparison Project Phase 5) for coupled model experiments 
(see: http://cmip-pcmdi.llnl.gov/cmip5/docs/Taylor_CMIP5_ 
design.pdf). 

As in PlioMIP Experiment 1, the concentration of CO 2 
in the atmosphere was set to 405 ppmv. In the absence of 
any adequate proxy data, all other trace gases and aerosols 
were specified to be consistent with the individual group’s 
pre-industrial control experiments. 

When trying to specify CO 2 as far back in time as 
the Pliocene epoch uncertainty is inevitable. Evidence for 
Pliocene CO 2 comes from a number of sources: (1) the stom- 
atal density of fossil leaves (Kiirschner et al., 1996), (2) car- 
bon isotope analyses (e.g. Raymo et al., 1996), (3) alkenone- 
based estimates (Pagani et al., 2010; Seki et al., 2010) and (4) 
boron isotope analyses (e.g. Seki et al., 2010). The values of 
CO 2 from each of these proxies differ. However, within error 
they overlap. The stomatal density records support a CO 2 of 
350 to 380 ppmv. The average of the Raymo carbon isotope 
analyses is similar to the stomatal-based estimates but peaks 
above that value (beyond 425 ppmv) occur and are entirely 
plausible. The Pagani et al. study looked at reconstructed 


CO 2 from a number of different marine records and the work 
shows clearly that in three of the six marine records a CO 2 
value of 400 or 405 is perfectly reasonable and the CO 2 range 
stated by Pagani et al. (2010) is 365 to 415 ppmv. In the 
Seki et al. (2010) study the alkenone-based CO 2 record is 
consistent with a value around 400 ppmv. The highest reso- 
lution section of the Seki et al. (2010) boron isotope record 
shows a decline in CO 2 from 400 to 280 after 3 million yr. 
Pliocene CO 2 is an important ongoing area of research with 
new records coming online in the next few years. Beyond Ex- 
periment 2, CO 2 is an obvious choice for sensitivity tests as 
part of PlioMIP and is one of the many justifications for util- 
ising coupled ocean-atmosphere models within the project. 

The solar constant and orbital configuration was specified 
as the same as each participating group’s pre-industrial con- 
trol run. The PRISM data set of mid-Pliocene boundary con- 
ditions (Dowsett, 2007) represents an average of the warm 
intervals during the time slab (3.29 to 2.97 million yr) rather 
than a discrete time slice, making it challenging to prescribe 
an orbital configuration which is representative of the entire 
~300 000 yr interval. Furthermore, it is difficult to provide 
an average insolation forcing at the top of the atmosphere in 
some climate models, with some models requiring specific 
values for eccentricity, obliquity and precession. Therefore, 
PlioMIP decided to specify a modern orbital configuration, 
even though available astronomical solutions (e.g. Laskar et 
al., 2004) indicate that this may not provide the most rep- 
resentative mean orbital forcing for the mPWP (see Hay- 
wood et al., 2010). In the future this uncertainty can be ex- 
plored through sensitivity experiments in which the orbital 
configuration is changed and provides a further justification 
of the implementation of coupled-ocean atmosphere models 
in PlioMIP. 

2.2 Adoption/availability of a “Preferred” and 
“Alternate” experimental design 

Two boundary condition data packages are available - “pre- 
ferred’’ and “alternate”. Versions 1.0 of both data pack- 
ages are provided on the PlioMIP website (http://geology.er. 
usgs.gov/eespteam/prism/prism_pliomip_data.html), and are 
provided as supplementary information to this paper. The 
preferred data package requires the ability to change the 
model’s land/sea mask to a mid-Pliocene configuration. The 
alternate data package, with a modern land/sea configura- 
tion, is provided in order to maximise the potential num- 
ber of participating groups in PlioMIP, since it is difficult in 
some climate models to successfully alter the land/sea mask. 
Groups that are not able to change their land/sea mask were 
asked to use their own modern land/sea mask. However, a 
PRISM3D/PlioMIP modern land/sea mask is provided in the 
alternate package to help guide the implementation of mid- 
Pliocene topography and vegetation etc. into different cli- 
mate models. 
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Table 1 . Experimental design - PlioMIP Experiment 2. 


Model Coupling 

Atmosphere-Ocean 

Integration Length 

500 yr (or as long as possible) 

Oceans 

Ocean Mode Deep Ocean Input 

Dynamic - initialized with PRISM3 ocean Global dot v2.0* if possible or 

temperatures if possible or Pre-Ind Control same as Pre-Ind Control 

Preferred Boundary Conditions 

Land/Sea Mask Topography Ice Sheets Vegetation 

PRISM3D PRISM3D PRISM3D PRISM3D or 

( land fraction vl.l) (topocvl.l*) (biome veg vl.3 or (biome _veg_vl. 3 or 

mbiome-veg-vl.3 ) mbiome-veg-vl.3) 

Alternate Boundary Conditions 

Land/Sea Mask Topography Ice Sheets Vegetation 

Local modern land/ PRISM3D PRISM3D PRISM3D 

sea mask ( topo_vl.4 *) (biome -Veg-V 1.2 or (biome _veg_v 1.2 or 

mbiome veg vl.l) mbiome veg vl.2) 

Greenhouse Gases 

C0 2 N 2 0 CH 4 CFCs o 3 

405 ppm As Pre-Ind Control As Pre-Ind Control As Pre-Ind Control As Pre-Ind Control 

Solar Constant 

As Pre-Ind Control 

Aerosols 

As Pre-Ind Control 

Model Spin-up 

Documented by individual groups 


* Applied as an anomaly to control experiment data sets used by each participating group rather than as an absolute. 


Those groups who are able to adopt the preferred exper- 
imental design will have to adjust the bathymetry of the 
deglaciated West Antarctic region. To avoid numerical in- 
stabilities groups will specify a flat bathymetry with an aver- 
age depth of 500 m. This bathymetry will be graded into the 
modern bathymetry. For those models which are very sen- 
sitive to bathymetric alteration, groups are free to do what 
is necessary to enable a simulation to be successfully com- 
pleted. 

Given the absence of a global data set for mPWP 
bathymetry in all other regions ocean bathymetry will be 
specified as modern. However, there has been a discussion 
recently regarding the depth of potentially critical sills in 
the North Atlantic (e.g. Jones et al., 2002). Robinson et 
al. (201 1) have shown that the exact depth of the Greenland- 


Scotland Ridge specified in coupled climate models can 
have a significant effect on model predictions of North At- 
lantic Deepwater formation (NADW), thermohaline circula- 
tion (THC) and increased SSTs in the North Atlantic and 
Arctic. Investigation of such uncertainties is beyond the 
scope of the first phase of PlioMIP but could become a fo- 
cus for sensitivity tests in the future. 

2.3 Implementation of ocean temperatures and 
topography as an anomaly 

To ensure that the climate anomalies (mid-Pliocene minus 
present day) from all PlioMIP climate models are directly 
comparable, i.e. that they reflect differences in the models 
themselves rather than the differences of modern boundary 
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conditions, it was decided to implement both the Pliocene 
topography and initial sea surface temperatures (SST) (from 
Experiment 1) and SST and deep ocean temperatures (from 
Experiment 2) as an anomaly to the standard modern ocean 
temperature and topographic data set used by each modelling 
group’s individual model. For Experiment 2 to create the 
initial Pliocene ocean temperature and topography, the dif- 
ference between the PRISM Pliocene and PRISM Modern 
ocean temperatures and topography will be calculated and 
added to the modern ocean temperature and topographic data 
sets each participating modelling group employs. 

In other words: 

Topo_Plio= (Topo_Plio_PRISM3D (1) 

— Topo_ModernJ :, RISM3D) 
+Topo_Moderniocal 
and 

OceanT_Plio= (OceanT_Plio_PRISM3D (2) 

— OceanT_Modern_PRISM3D) 
+OceanT_Modern_Local 

However, when using such a method, a potential mis- 
match between mid-Pliocene and modern topography land- 
sea masks is possible. This will be overcome by using abso- 
lute Pliocene topography and ocean temperatures in regions 
where no modern data are given (such as for the Pliocene 
topography in the Hudson Bay region). Modern SSTs are 
projected on the same Pliocene grids {preferred and alter- 
nate) to make anomalies easier to generate. There may be 
mismatches (for example in the West Antarctic region) be- 
tween the Pliocene deep ocean temperature data, where it is 
provided, and the Pliocene land/sea mask. Where this is the 
case. Global jdot_v2.0 deep ocean temperatures should be ex- 
trapolated horizontally into regions with no data coverage, 
therefore maintaining the Global jdot v2.0 vertical tempera- 
ture profile. Groups unable to alter the initial ocean tempera- 
ture state should begin the simulation with a modern control 
state and document the spin-up of the simulation. Salinity 
should be derived from Levitus and Boyer (1994) or an ex- 
isting modern (control) simulation. The starting atmospheric 
conditions for Experiment 2 should be derived from the end 
of Experiment 1 if possible. 

3 Description of Boundary Conditions (PRISM3D) 

A full description of the mPWP land-sea mask (including 
the nature of ocean gateways) and topography (outside of ice 
sheet regions), ice sheet height and extent, SST (see Fig. 1), 
sea-ice extent, vegetation type and distribution, soils, lakes 
and river routing is provided in Haywood et al. (2010; this 
volume). Here it is only necessary to describe the construc- 
tion and nature of the three-dimensional data set for ocean 
temperatures that groups have the option of using to initialise 
their ocean models for PlioMIP Experiment 2. 


3.1 The PRISM3D data set of ocean temperatures 

The PRISM3D deep ocean temperature reconstruction 
(Fig. 2) is presented at a 4° latitude by 5° longitude reso- 
lution with 33 depth layers, based upon 27 localities, un- 
evenly distributed among the ocean basins. While not op- 
timal for generating a global reconstruction, this represents 
possibly the largest number of temperature estimates for 
any deep-water reconstruction from any time interval. Spe- 
cific steps in the reconstruction methodology are listed in 
Dowsett et al. (2009) Supplement A (http://www.clim-past. 
net/5/769/2009/cp-5-769-2009-supplement.pdf). Because 
the PRISM3D reconstruction is designed for coupled ocean- 
atmosphere general circulation models, it represents a recon- 
struction of a prescribed day, arbitrarily chosen to be 1 De- 
cember (this does not imply that climate models must be in- 
tegrated from this date). The PRISM3D November and De- 
cember monthly SST reconstructions were averaged to ap- 
proximate the SST for mid-Pliocene 1 December. A surface- 
temperature anomaly was created by subtracting the modern 
1 December SST field (Reynolds and Smith, 1995) from the 
Pliocene 1 December data. The surface temperature anomaly 
was then added to the 0 m layer of Levitus and Boyer (1994) 
(converted to a 4° x 5° resolution) to create the PRISM3D 
0 m reconstruction. Since no data points fall between 0 m and 
1 100 m in the deep ocean temperature data set, PRISM chose 
to use a mathematical function that decreases the weight of 
the surface anomaly with depth down to 1400 m (see Dowsett 
et al., 2009). Between 900 m and 1400 m that anomaly was 
further modified based upon data from Southern Ocean sites 
to accomplish an adjustment or vertical expansion of palaeo 
Antarctic Intermediate Water (AAIW). 

In the Atlantic sector of the Southern Ocean, data suggest 
warmer palaeo NADW expansion in the region relative to 
modern day (Dowsett et al., 2009). Warm anomalies for the 
mid-Pliocene at all sites are in keeping with the hypothesized 
warmer and stronger flux of palaeo NADW and diminished 
(colder) palaeo Antarctic Bottom Water production relative 
to today. In the western Pacific, sites that monitor Pacific 
Deep Water today show small positive anomalies for the mid- 
Pliocene. PRISM interprets these data to indicate the overall 
warmer conditions of the water masses that mixed to form 
palaeo PDW. In the eastern Pacific, temperature data can be 
explained by a vertical displacement of AAIW concomitant 
with overall warming of PDW. 

A complete discussion of the rationale and methodology 
used for the PRISM Deep Ocean Temperature reconstruction 
can be found in Dowsett et al. (2009). 


4 Data management and planned analyses 

PlioMIP has adopted the established variables list out- 
lined by the second phase of the PMIP (Braconnot 
et al., 2007a, b). Model outputs will be submitted 
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February SST Anomaly (PRISM3 - Modern) August SST Anomaly (PRISM3 - Modern) 



PRISM3D Sea Ice Extent (February) PRISM3D Sea Ice Extent (August) 



Fig. 1. PRISM3D SST anomaly for February (top left) and August (top right). PRISM3D sea-ice extent for February (bottom left) and 
August (bottom right). 



Fig. 2. Longitudinal profiles of ocean temperature from transects at (a) 15° W, (b) 45° W, (c) 90° W and (d) 165° E. All temperatures are 
shown in °C. Contour interval is 2 °C. Black contour lines show modern temperature overlaid on coloured regions showing the mid-Pliocene 
reconstruction. The change in temperature can be surmised by comparing the colour contours to the black overlaid contour lines. For 
example, in panel (d), Site 806 is slightly cooler than 2 °C in the modem ocean but was slightly warmer than 2 °C during the mid-Pliocene 
(modified from Dowsett et al., 2009). 


www. geosci-model-dev. net/4/57 1/2011/ 


Geosci. Model Dev., 4, 571-577, 2011 







576 


A. M. Haywood et al.: Pliocene Model Intercomparison Project (PlioMIP) 


and stored within the PMIP2 database. Specifically, 
for PlioMIP Experiment 1, this refers to PMIP2 rec- 
ommended outputs for the atmosphere (outlined on 
the PMIP2 website http://pmip2.lsce. ipsl.fr/>Experimental 
Design> Variables > Atmosphere). PMIP/PlioMIP requires 
participants to prepare their data files so that they meet the 
following constraints (regardless of the way their models pro- 
duce and store their results). 

- The data files have to be in the (now widely used) 
netCDF binary file format and conform to the CF (Cli- 
mate and Forecast) metadata convention (outlined on 
the website http://cf-pcmdi.llnl.gov/). 

- There must be only one output variable per file. 

- For the data that are a function of longitude and lati- 
tude, only regular grids (grids representable as a Carte- 
sian product of longitude and latitude axes) are allowed. 

- The file names have to follow the PMIP2 file name con- 
vention and be unique. 

Participants are encouraged to create the files for submis- 
sion to the database using the CMOR (Climate Model Output 
Rewriter) library. This library has been specially developed 
to help meet the requirements of the Model Intercompari- 
son Projects. Details of the CMOR library are provided on 
the PMIP2 website (http://pmip2.lsce. ipsl.fr/>Experimental 
design> Output format>CMOR library). Proposals for 
model analyses using PlioMIP Experiment 1 data can be 
made using the established protocols outlined on the PlioMIP 
website. 

5 Conclusions 

This paper provides a model intercomparison project de- 
scription for the Pliocene Model Intercomparison Project 
(PlioMIP) and documents in detail the experimental design. 
Specifically, this paper describes the experimental design and 
boundary conditions utilised for Experiment 2 of PlioMIP, 
following a companion paper for Experiment 1. Experi- 
ment 2 will utilise coupled ocean-atmosphere models and 
will enable the sensitivity of the simulated coupled ocean- 
atmosphere system in a palaeoclimate context to be explored. 
It also provides the necessary foundation for further sensitiv- 
ity tests as part of PlioMIP that will explore uncertainties in- 
troduced in the climate modelling from specified trace gases 
and orbital configurations etc. 

Supplementary material related to this 
article is available online at: 
http ://www.geosci- model- dev.net/4/57 1/2011/ 
gmd-4-571-2011-supplement.zip. 
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